Eecient Support of Parallel Sparse Computation for Array Intrinsic Functions of Fortran 90 *

Fortran 90 provides a rich set of array intrinsic functions. Each of these array intrinsic functions operates on the elements of multi-dimensional array objects concurrently. They provide a rich source of parallelism and play an increasingly important role in automatic support of data parallel programming. However, there is no such support if these intrinsic functions are applied to sparse data sets. In this paper, we address this open gap by presenting an eecient library for parallel sparse computations with Fortran 90 array intrinsic operations. Our method provides both compression schemes and distribution schemes on distributed memory environments applicable to higher-dimension sparse arrays. This way, programmers need not worry about these low-level details. Sparse programs can be expressed concisely using array expressions, and parallelized with the help of our library. Our current testbed is an IBM SP2 workstation cluster. Preliminary experimental results show that our approach is promising in speeding up sparse matrix computations on both sequential and distributed memory environments if the computations are expressed with Fortran 90 array expressions.